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The forward scattering pattern from a transversely illuminated 
polymer-coated optical fiber is modeled with a ray- tracing program. 
For polymers with a refractive index higher than that of silica, the 
scattering pattern contains caustic surfaces. If the coating is eccentric 
in a direction perpendicular to the incident rays, we show that the 
positions of the caustic surfaces (6a) vary monotonically with eccen- 
tricity at the rate of 0.85° /fim for small eccentricities for a coating 
with a 1.539 refractive index and 59.3-iim thickness on a 110-fim 
diameter silica fiber. Another identifiable feature of the pattern (6 b) 
has a predicted small-signal sensitivity of 1.1° /fim. The information 
obtained from the pattern can be used to insure well-centered coatings 
on optical fibers. A simple interpretation of the scattering pattern 
provides an error signal to guide compensatory coating cup align- 
ment. The results are reported of well-centered runs obtained on two 
separate drawing lines by this means. 

I. INTRODUCTION 

The goal of placing an axially symmetric coating around an optical 
fiber has been experimentally approached through various tech- 
niques 1 " 4 and theoretically analyzed with various models. 4-7 There is 
ample reason for this interest. Good centering has been shown to 
increase fiber strength. 8 It also eases the fabrication of connector 
arrays for splicing and provides maximum abrasion protection for a 
given coating. Furthermore, an asymmetric coating may add to fiber 
loss. 29 

For a coating that is fairly transparent and has a refractive index 
more than that of the fiber cladding, a light beam incident essentially 
transverse to the optic axis has been found to have a notably charac- 
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teristic forward scattering pattern. Its shape and position are easily 
monitored and provide unambiguous information about both coating 
thickness and eccentricity for any degree of eccentricity. The back- 
scatter pattern has been previously analyzed and used as an inline 
centering monitor. 3,7 

In this paper, we describe the pattern, recount our experiences with 
it when used as a coating monitor, and present a ray trace model of 
the forward scattering. We use this model, together with a more limited 
closed form expression, to obtain a better understanding of the sensi- 
tivity and operating range of this monitoring technique. The pertinent 
parameters are coating refractive index, coating thickness, and eccen- 
tricity. 

II. THEORY 

Under appropriate values of coating refractive index, thickness, and 
eccentricity, an optical fiber illuminated as in Fig. 1 generates the 
scattering pattern shown in Figs. 2 and 3. We wish to predict the 
scattering angles A and 6 B associated with the A and B boundaries 



COATED FIBER 
SCREEN 

SCATTERING 
f PATTERN 

I 

\ 




Fig. 1— Setup for generating the forward scattering pattern of a coated fiber off-line. 
A cylinder lens is used here to vertically expand the incident laser beam. This allows the 
coating geometry of a stationary fiber to be directly compared from point to point over 
the height of the pattern. Such vertical resolution may prove to be a useful tool for 
diagnosing coating defects (Ref. 3). 
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Fig. 2— Plot of the envelope of the characteristic scattering pattern intensities for (a) 
centered and (b) eccentric coatings (interference fringes ignored— see Fig. 3). There is 
a central bright zone between the points labeled A and another, but much less intense, 
bright zone beyond each B point. 

between light and dark zones. Of special interest are the locations of 
the A boundaries. Their clear delineations enhance their potential as 
reference points for automated derivation of error signals. 

Since our interest is primarily to predict the angles A and 0b, a 
prediction of the power envelope of the scattering pattern (as opposed 
to individual fringes) should suffice. This approach allows a variety of 
simplifications and the model thereby obtained can be characterized 
by the following list of assumptions: 

1. Ignore phase and polarization of scattered radiation. 

2. Assume both fiber and coating to be cylindrical, though not 
necessarily concentric. 

3. Allow eccentricities only transverse to the incident rays. 

4. Ignore Fresnel reflections. 

5. Assume lossless dielectrics. 

6. Assume n = 1.457 at 633 nm throughout cladding and core (ignore 
core profile). 

7. Assume incident rays are all parallel and of uniform intensity. 

As we shall see later, the model's predictions of A and 0n agree 
closely with the experiment. However, to appreciate its realm of 
applicability, some point-by-point justifying comments are helpful. 
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Fig. 3— Photograph of a scattering pattern from a cured, n = 1.539 coated fiber. 
Interference fringes can be seen (not included in Fig. 2). A and B indicate bright-dark 
boundaries as in Fig. 2. C marks that part of the incident beam that was not scattered 
by the fiber but continued straight on to the screen. D indicates an out-of-focus image 
of part of the fiber support and is unrelated to the pattern seen on the screen. 

(i) The inclusion of phase and polarization information in the 
scattered radiation would enable us to predict the structure of 
the interference fringes as was done by Watkins. 10 But our 
current concern for merely the power envelope of the pattern 
led us to ignore these considerations. (If an automated scatter- 
ing monitor is unavoidably affected by the interference pattern, 
then a refined, phase-sensitive model would be required for 
design guidance.) 

(ii) In the great majority of the coated fibers, both coating and 
fiber are nearly circular cylinders. 

(Hi) While allowing only transverse eccentricity does restrict the 
applicability of the results, this is not troublesome from an 
instrumentation point of view. The inline monitoring equip- 
ment is always an orthogonal cross-beam set-up (see Section 
V), so that aligning for symmetry on one axis assures only 
transverse eccentricity on the other. 

(iu) Except for relatively few rays (near total internal reflection 
conditions), the Fresnel reflections carry too little power to 
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materially alter the high irradiance gradient observed at the 
boundaries, particularly the A-boundaries. 
(v) By reasoning similar to (iv), a small transmission loss through 

one pass of the coating can be ignored. 
(vi) Rays passing through the fiber core all lie within the bright 
central region of the scattering pattern and do not affect the 
contrast at A or B. 
(uii) The incident HeNe laser beam provides a reasonably plane 
and uniform intensity source over the <0.3-mm outside diam- 
eter of the coating. 
The model is embodied in a ray-tracing Fortran computer program, 
vscate. Each ray is identified by its normalized height H, the height 
of the ray above the coating's center expressed as a fraction of the 
coating radius r, (see Fig. 4). After tracing the ray, we find it exits the 
coating at some angle, 6, relative to the incident direction. From the 
discrete pairs of values (0, H), we can approximate dd(H)/dH. The 
uniform incident radiation assumption then allows us to state 



P(d)c*l\dd(H)/dH\ 

H 



(1) 




Fig. 4— The three characteristic types of ray paths: P — penetrates the fiber, tir— 
undergoes a total internal reflection at the coating-fiber interface, and CO — passes 
through the coating only. Fiber radius is 55 ftm, coating radius is 114.3 fim, and an 
eccentricity (E) of 10 /im exists. Refractive indices are no = 1.0003, rt\ = 1.539, and n-i 
= 1.457. Sign conventions for E, H, and scattering angles (0) are indicated by the arrows. 
H- values are 0.40, 0.57, and 0.65 for the P, tir, and co rays, respectively. 
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where P(0) is the power scattered per unit angle at 6. The summation 
is required because more than one H may correspond to a particular 
6. For tracing rays incident below the coating's center (H < 0), we 
appeal to the symmetry 

0(H, E) = - 8(-H, -E), (2) 

where E is the signed eccentricity of the coating. The divergence of 
eq. (1) when dO/dH = is a consequence of ignoring diffraction effects 
and is related to assumption 1 of the model. The integrated scattered 
power does not diverge, whether or not diffraction is included. 

In Fig. 4, we have a cross section of a coated fiber with three rays 
traversing it, showing the sign conventions for H, E, and 0, and the 
three characteristic paths a ray may take once it enters the coating: 
(i) passes through the coating and penetrates (P) the fiber, (ii) expe- 
riences a total internal reflection (tir) at the coating-fiber interface, 
and (Hi) traverses the coating only (co) and does not intercept the 
fiber. (Actually, a fourth possible path exists for sufficiently eccentric 
coatings in which a tir occurs as the ray tries to exit the coating. This 
case is of significance only for large eccentricities as will be seen in Fig. 
7.) The three rays illustrated were calculated by the vscate ray trace 
program with no = 1.0003, /ii = 1.539, and n 2 = 1.457, appropriate for 
a coating with a refractive index of 1.539 on a solid 100-percent silica 
fiber, surrounded by air. 

Ray paths are uniquely determined by the if-value, refractive indi- 
ces, and cross-section geometry. An equation representing each 
straight-line segment is obtained by appropriate application of either 
Snell's Law or a tir at each preceding boundary. A simultaneous 
solution of this equation and the next boundary locus provides the 
next boundary intersection, whereupon the process repeats. Finally, as 
the ray exits the coating, its overall deviation from its original 
direction is related to the original H. 

For the concentric coating case, 6(H) can be obtained in closed form 
by application of Bouguer's formula" without a ray trace. Details are 
given in Appendix A. A Fortran program, scatcf, was written to 
evaluate the scattering pattern predicted by the closed form expres- 
sions. Although applicable only to the concentric case, it is valuable 
when compared to the more general ray trace approach as a check for 
mutual consistency between the models. Excellent agreement was 
found. For example, Figs. 8 and 9 can be obtained identically from 
either vscate or scatcf. Also, the relatively simple form of scatcf 
shows better the operating limits of the monitoring system in terms 
of coating thickness and refractive index. For example, if n 2 r 2 //iori > 
1, then all arguments in eq. (lie) are less than unity for all values of if. 
We interpret such a condition to mean that eq. (lie) is valid for all 
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rays; that is, all rays penetrate the glass (P-rays). Similarly, if n 2 r 2 / 
non < 1 and /iir 2 /Wi > 1, then eq. (lib) is valid for some H; that is, 
in addition to fiber-penetrating rays, some rays experience a tir at the 
plastic-fiber interface (TiR-rays). But there are still no rays of the 
coating-only type (co-rays). Finally, if n^/rton < 1, then all types of 
rays are present and a full scattering pattern forms. Table I gives the 
conditions for the existence of the various rays and the values of H 
corresponding to each type. 

III. PREDICTIONS 

Plots of 6{H ) are perhaps the clearest way to appreciate the origin 
of the light and dark zones. The solid line in Fig. 5 plots 0(H) for a 
fiber with a cured concentric n = 1.539 coating with n = 114.3 /mi and 
r 2 = 55.0 /tm. The assumed refractive index values are the same as 
those assumed in drawing the Fig. 4 rays. Upon comparing this result 
to the scattering of a homogeneous cylinder of n = 1.539 coating of 
114.3-/mi radius (long-dashed line), we note the following: (i) The two 
plots differ most dramatically at H = /i 2 r 2 //ion (~0.70). It is at this H- 
value where the silica cylinder is acting most strongly as a negative 
lens; ray refractions or tirs at the fiber-coating boundary are most 
strongly counteracting the ray refractions at the coating-air boundary. 
Also, at this H, the transition from P to tir rays occurs, (ii) The zero 
value of dd/dH at H -0.62 implies the formulation of a caustic surface 8 
at the corresponding angle, 18.45 degrees. The bright caustic surface 
is represented by the peaks in Fig. 2. (Hi) At slightly greater scattering 
angles, light comes from rays originating at H ~0.73, where dd/dH is 
relatively large, implying a low P(0). We therefore predict a high 
contrast between a bright inner zone (ending at 18.45 degrees) and the 
angles immediately beyond, just as observed at boundary A. Starting 
at if = nir 2 /nori (0.74), corresponding to at 38°, the rays miss the 
fiber and traverse only the coating. As it should, the plot thereafter 
coincides with the coating-only case. Also, the sudden decrease in dd/ 
dH occurring at 38 degrees predicts boundary B, the start of the outer 
bright zone. 

We see in Fig. 5 how a change in coating refractive index to 1.49 

Table I — Conditions for the presence of the various 
ray types 

Conditions on n,,, n\, n>, n, r-i 
Types of Rays Present When E = 

P n*n/nori > 1 

Pand tir nir./nur, < 1 and n l r>/n { ,r l > 1 

P, TIR, and CO Rira/flori < 1 
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Fig. 5— Plots offfvsH for three profiles. In all cases, no = 1.0003 and n 2 = 1.457. For 
the longer dashes, r 2 /r, = (no fiber, merely a cylinder of coating) and n\ = 1.539. For 
the solid line, r 2 /n = 0.4812 and n t = 1.539. For the shorter dashes, r 2 /r t - 0.4812 and 
n, = 1.490. 
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changes 0(H) (short dashed lines). The key difference is that the broad 
peak has narrowed, indicating a decrease in illuminance at the edge of 
the bright zone. The co rays do not superimpose on the n = 1.539 
homogeneous cylinder, since the refractive power is less for the lower 
index. 

Figure 6 indicates how an eccentricity changes the pattern. With 
dimensions unaltered from Fig. 5, an eccentricity of 10 fan transverse 
to the incident rays is introduced. We see that the scattering pattern 
is now asymmetric about H = and its orientation relative to the 
eccentricity is shown in the drawing. 

By evaluating a variety of eccentricities, the sensitivity of the 
scattering pattern to the transverse offset can be ascertained. 6a and 
0b are plotted against eccentricity in Fig. 7 for left and right sides of 
the scattering pattern. We see that, for rii = 1.539 and r 2 /ri = 0.4812, 
one set of B A and 0b boundaries disappears at E = ±22 A and ±26.2 nm, 
respectively, while the other set moves toward the center and remains 
visible even when the fiber reaches the maximum possible eccentricity. 
The $ A disappearance is a result of the previously mentioned tirs, 
which can occur as the ray attempts to exit the coating while the Ob 
limit, in this case >90°, is a manifestation of the maximum refracting 




t "fiber ""*J 



Fig. 6 — Plot of vs H for the case of Fig. 4. Spatial relationships among E, H, and 8 
are depicted in the inset. The plotted points are computed. 
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Fig. 7 — Plots of 6a and do on both sides of the scattering pattern as they would appear 
to an observer as the Fig. 4 fiber coating became eccentric to the left (up) and to the 
right (down). 

power of the coating. Error signals for an automated system are in 
principle obtainable from the non-disappearing 0a for all eccentricities. 
The computed slopes of the Fig. 7 curves at E = show a small- 
eccentricity sensitivity of 0.85 degrees/ /im for the two A 's and 1.1 
degrees/jtim for the two 0b's. 

The angle separating B from A varies with n\ as other parameters 
are held constant. Figure 8 illustrates this with plots of 0a and Ob as 
functions of rti. 0a and B also vary significantly with the ratio r 2 /r\ , as 
shown in Fig. 9. This dependency can conceivably be used for moni- 
toring coating thickness inline, given the fiber diameter. Instrumen- 
tation for this would be more complex than that used for merely 
monitoring eccentricity. Notice that the 0a caustic disappears when 
r 2 /ri is below a certain value (-0.10 in Fig. 9). When r 2 /ri is sufficiently 
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Fig. 8—6 A (n,) and H (n,) for the case n„ = 1.0003, n 2 = 1.4570, r 2 /r, = 0.4812, and 
£ = 0. 

small, the negative lensing at the coating-fiber interface is more 
powerful than the positive lensing of the air-coating interface, even for 
paraxial rays (\H\ <sc 1). That is, the condition that the converging 
wavefront formed by the air-coating interface be refracted to a plane 
wavefront at the coating-fiber interface leads to the lower 0a disap- 
pearance condition 



ri _ no(rii — 712) 
n n.2(ni — no) 



(3) 



derived in Appendix B. 
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Fig. 9— 6 A (r>/n) and 6 u (r 2 /n) for the case n„ = 1.0003, n, = 1.539, n 2 = 1.457, and 
E = 0. In general, 6a exists (a caustic forms) when n n (n\ — n-i) / n-i(n\ — no) < r-i/r, < ru,/ 
n.2 and 6b exists when ra/fi < ito/nj. 



This evaluates to r 2 /ri = 0.1045 for an n = 1.539 coating which is 
the lower 0a cutoff point, as seen in Fig. 9. In actual fibers, the presence 
of a higher-index core will tend to slightly decrease the cutoff ratio 
below that predicted by eq. (3). 

IV. COMPARISON WITH EXPERIMENT 

The scattering pattern from the well-centered n = 1.539 coating 
shown in Fig. 10 was examined in detail. The fiber was rotated through 
360 degrees of azimuth in 10-degree increments, and measurements of 
the positions of the respective A and B boundaries were made on a 
screen 47 mm behind the fiber. The wavy nature of the pattern on the 
screen, seen in Fig. 3, gave each measurement a range of values 
indicated by the width of the lines in Fig. 11. A statistical summary is 
given in Table II. 

With values of no = 1.0003, n 2 = 1.539, n 2 = 1.457, and r 2 //*i = 0.553, 
we compared this observed pattern to the theoretically predicted one. 



324 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1980 




CENTERING OPTICAL FIBER COATINGS 325 



SCREEN POSITION 



no 



B LEFT 

50 40 



A 
30 20 




20 
40 

60 I- 
80 - 



• mm . 




A 
20 30 



RIGHT 
40 



B 

50 



i 100 





120 


< 









IT 


140 


N 




> 




< 


160 


ir 




t- 




:d 


im 






< 








z 
o 


200 






< 




5 


220 



O 240 



Z 260 



280 
300 
320 
340 h 



360 



J «L 



51.9 46.8 40.4 32.6 23.1 12.0 
9 B A — 

LEFT 



12.0 23.1 32.6 40.4 46.8 51.9 



Fig. 11 — Values for A and B on both sides of the scattering pattern as the Borden- 
coated fiber of Fig. 10a was rotated in 10-degree increments. The width of lines 
represents the position uncertainty due to readability and pattern waviness (see Fig. 3). 
The screen was 47 mm behind the fiber. Statistics are given in Table II. 



They were mutually consistent: The average measured value of the 
A was 25.7° ± 2.7° versus 25.4° ± 1.3° predicted from the ray trace. 
For Ob, we measured 45.7° ± 3.2° versus 49.4° ± 2.2° predicted. The 
estimated experimental errors reflect both the distribution in measured 
values from different azimuths and uncertainty in each single mea- 
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surement. Theoretical uncertainties arose from inexact knowledge of 
ri and r2 as obtained from cross sections. 

V. USE AS A COATING MONITOR 

The pattern's high sensitivity coupled with its easy identifiability 
fill the requirements of a coating geometry monitor very well. A 
skewed pattern is interpreted as an error signal. Then, by monitoring 
scatter patterns from two perpendicular directions, one can guide the 
coating applicator into proper, centered alignment. Continued moni- 
toring and feedback would compensate for subsequent drifts in the 
fiber path. 

Figure 12 shows a crossed beam setup used to obtain the centered 
n = 1.539 coating over a silica (T08) filament shown in Fig. 10a. Figure 



Table II — Comparison of experimental results with 
theoretical predictions 



A (mm) 



fl(mm) 



Right side 
Left side 
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Theoretical 
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Fig. 12— Plan view schematic of the setup used inline on the Bell Laboratories, 
Atlanta, drawing machine to monitor coating geometry. The beams are slightly offset 
vertically to avoid ambiguities on the screen. 
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10b shows the original eccentric condition before alignment. Vertical 
beam spread was omitted here and the scattering pattern was a bright 
modulated horizontal line. An n = 1.49 coating was also successfully 
centered, but the lower refractive index of the coating decreased 
pattern visibility as compared to the higher-index coating. 

A similar crossed-beam monitoring setup was then used on another 
drawing line to monitor and maintain centering on a 6.6-km length of 
optical waveguide. Table III shows how eccentricity varied with posi- 
tion as measured on cross sections taken later. Overall, we see excellent 
centricity, but at the 4.4-km position an unexplainable gross eccentric- 
ity occurred. It was traced forward and back and found to have started 
about 250 m before and ended 100 m after the 4.4-km point. Nothing 
in the ray trace theory suggests any ambiguity in the pattern at large 
eccentricities. 

In the past, unpublished data had been taken by the author on the 
results of using alignment techniques which did not utilize light scat- 
tering. The best set of data showed a median eccentricity of 12.9 pan 
with other medians from about 15 /tm to 25 jum. Also, these other 
techniques could not detect ongoing changes in alignment during the 
drawing process. 

VI. SUMMARY 

We have presented a ray trace model which indicates that the 
forward scattering pattern from a transversely illuminated optical fiber 
provides unambiguous data about the eccentricity of the coating under 
certain conditions of fiber size, coating thickness, and coating refractive 
index. 

Furthermore, the data can be easily interpreted to provide an 
eccentricity error signal. The ray trace model has been shown to agree 
well with a closed-form expression of more limited applicability. It also 

Table III — Coating geometry of a 6.6-km optical fiber where forward 

scattering monitoring was used to control coater alignment. 

Monitoring and feedback started after the 0-km point. The excursion 

at 4.4 km remains unexplained. Uncertainty of all cross-section 

measurements is within ±2 /xm. 



Distance from 


Coating 


Glass 


Coating 


Start 


Eccentricity 


Diameter 


Diameter 


(km) 


(/mi) 


(/mi) 


(/mi) 


0.0 


46.3 


110.6 


245.6 


1.1 


5.0 


109.3 


252.4 


2.2 


7.7 


108.7 


260.7 


3.3 


2.0 


109.8 


265.8 


4.4 


59.5 


110.6 


256.8 


5.5 


8.1 


110.5 


268.2 


6.6 


9.7 


108.7 


280.9 
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matches scattering data from experiments. Finally, when the forward 
scattering pattern was used to monitor coating geometry inline with 
the drawing and coating process, the information it provided allowed 
us to apply well-centered coatings. 
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APPENDIX A 

The Concentric Case 

An optical fiber and its concentric coating have a refractive index 
dependent solely on the radial distance from their common center, 
that is, 

n = n(r). (4) 

In such a case, Bouguer's formula 11 applies: For a ray lying in a plane 
perpendicular to the fiber's axis, the product of the local refractive 
index and the perpendicular distance to the fiber's center is constant 
along the ray's path, that is, 

n(r) d(r) = C, (5) 

where C is the constant and d(r) is the perpendicular distance from 
the tangent to the ray at r to the center (r = 0) (see Fig. 11). From this 
relationship and symmetry arguments, we derive the dependence of 
the scattering angle, 6, on the normalized height, H. 

First we note, in Fig. 13, the point A where an arbitrary ray BAC, 
with normalized height H, has traveled half its total distance through 
a cylinder where ti = n(r). We observe a mirror image symmetry for 
the ray path with respect to the straight line DAF containing the 
halfway point A and the fiber center F. Since half the refraction occurs 
on each half of the ray path, then the ray path has turned through an 
angle 0/2 by the halfway point. Furthermore, this halfway point is the 
point of minimum r for the ray (assuming some other minimum point 
leads to either a symmetry violation and/or a contradiction with 
Bouguer's formula). 

In terms of Fig. 13 nomenclature, we have (Ref. 11, p. 123, eq. (10)); 

f rmin Cdr 
d <t> = , (6) 
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Fig. 13— Nomenclature for derivation of closed form scattering expression from 
Bouguer's formula (see Appendix A). 



Also, from Fig. 11, we have the geometric relations 

6(H) = -n-2$ (7) 

and 

cos(<fo - 0/2) = H. (8) 

For the ray in Figure 11, 

C-noflh. (9) 

Substituting eqs. (7), (8), and (9) into eq. (6) and integrating the lhs 
gives 



0(H) = -2 cos-'iH) + 2//rin<, 






nlH 2 r\ 



(10) 



For the two cases, (i) ray bypassing the fiber and traversing the coating 
only and (ii) ray penetrating the fiber, the second term on the right of 
eq. (10) evaluates to different forms (n(r) makes a step at the coating- 
fiber interface, r ™ Ma). For a third case, the ray experiences a total 
internal reflection (tir) at the coating-fiber interface, r = r2. There, 
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the turn in the path due to the tir must be subtracted from 0(H). 
After some algebra, 12 the resultant closed forms are 



2[sin- , if-sin- I -tf 



21 sin" 1 H - sin -1 — H + sin ' — '-'- H 



0(H) = < 



n ] 






21 sin" 1 H - sin" 1 — H + sin -1 "" - H 



n\ 



., Ho r x 
n x r 2 



(Ha) 
(lib) 
(lie) 



- sin H 

n 2 r> 

where (11a) applies to the coating-only rays (n\r 2 /noT\ < H < 1), (lib) 
applies to the tir rays (n-irz/riur) < H < /iir 2 //iori), and (lie) applies 
to the fiber penetration rays (0 < H < tiir-i/norx). 

These expressions were then incorporated into a Fortran program, 
scatcf, to determine 6(H) and d6/dH for the concentric coating case. 

APPENDIX B 

Derivation of Equation (3) 

We wish to find a relationship among the radii and refractive indices 
to generate the behavior of the incident wavefront J as shown in Fig. 
14: The incident flat wavefront J is refracted at the cylindrical air- 
coating boundary of radius r\. In the paraxial approximation, the 
refraction forms a cylindrical wavefront K of radius pi . It then proceeds 
a distance n — r 2 to L where it has a radius p 2 where 

P2 = pi - (n - r 2 ). (12) 



At this point it undergoes a second refraction at the coating-fiber 
interface at radius r 2 to again assume the flat shape shown at M. This 
wavefront behavior, where the paraxial refracting power of the inner 
cylinder just cancels the outer cylinder's is the limiting case for the 
disappearance of 0a (see Fig. 9). 




Fig. 14 — Nomenclature for derivation of lower cutoff value of r-i/r\ for existence of 
9a (see Appendix B). 
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The relationship among ?to, fii t n and p\ is easily derived and given 
in many texts on geometrical optics. It is a special case of the Gaussian 
formula for a single spherical surface 13 (or cylindrical surface for this 
two-dimensional case). Applied to this case we have 

pi -— i— . (13) 

1 - (no /ni) 

The flattening of the wavefront at the second surface is equivalent 
to a flat front entering from the right and transformed to a front 
converging to the left with radius p2. Therefore, again using the 
Gaussian formula, 

P2=- r"TT- (14) 

1 - (n 2 /ni) 

Combining eqs. (12) to (14) and eliminating p\ and p 2 , we arrive at 

r<i noini — m) ,,_. 

— =—, r. ( 15 > 

n 7i 2 (m - no) 

which is eq. (3) in the main text. 
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